Bang-Bang control of a qubit coupled to a quantum critical spin bath 
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We analytically and numerically study the effects of pulsed control on the decoherence of a qubit 
coupled to a quantum spin bath. When the environment is critical, decoherence is faster and we 
show that the control is relatively more effective. Two coupling models are investigated, namely a 
qubit coupled to a bath via a single link and a spin star model, yielding results that are similar and 
consistent. 
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I. INTRODUCTION 



Decoherence results from the unavoidable coupling be- 
tween any quantum system and its environment, and is 
responsible for the dynamical destruction of quantum su- 
perpositions, ft is detrimental for quantum information 
processing [T| since it leads to a loss of the quantum par- 
allelism that is implicit in the superposition principle. 
The possibility of preventing or avoiding decoherence is 
hence of significant importance for any technological use 
of quantum systems, aimed at processing, communicat- 
ing or storing information. To this end, one must un- 
derstand and model all of the relevant features charac- 
terizing the environment of the physical system to be 
protected. Understanding decoherence is also of funda- 
mental interest in its own right, since it is at the basis of 
the description of the quantum-classical transition Q ■ 

The study of open quantum systems has a long history, 
and many ingenious models have been proposed in order 
to describe the action of the environment in a quantum 
dynamical framework (see, e.g., Ref. [|[). Paradigmatic 
models represent the environment as a set of harmonic 
oscillators [1] or spins [f|. Recently there has been a re- 
newed interest in the analysis of decoherence induced b> 
such spin baths |6, 
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1231 [24| ; these are clearly relevant 
in a number of physically important situations, such as 
NMR [25j or spin qubits (26[, where loss of coherence 
is induced by the coupling to nuclear spins |27|. Sev- 
eral questions have been addressed so far and the picture 
that emerges is rather rich. A possible, monogamy-like, 
relation between the entanglement in the bath and de- 
coherence has been put forward in Ref. [ll[ and subse- 
quently analyzed in different papers. The signatures of 



criticality of the environment in decoherence have been 
discussed throu gh t he study of solvable one-dimensional 
model systems jl3l . EH EI|- A universal regime exists, 
in the strong coupling limit, in which the decay of the 
Loschmidt echo |28| does not depend on the system-bath 
coupling [Hi!]. 

Several different protocols have been designed to pro- 
tect quantum information. These include passive cor- 
rection techniques, in which quantum information is en- 
coded in such a way as to suppress the coupling with 
the environment [29|, |30( , and active approaches such as 
quantum error correction [l| and dynamical decoupling 
techniques [U [H, [H, (for an overview see, e.g., [HI], 
for historical references see [36j]). Dynamical decoupling 
strategies aim, by means of a dynamical control field, at 
averaging to zero the unwanted interaction with the envi- 
ronment. In its simplest version, which we consider here, 
the control field comprises a train of instantaneous pulses 
( "bang-bang" control) . While previous work on dynami- 
cal decoupling has made clear distinctions between differ- 
ent environments, in particular bosonic baths [31[ versus 
spin baths [13, HI], and fast versus slow 1// noise [39 ]. no 
attention has been paid so far to the impact a quantum 
critical environment might have on the efficacy of decou- 
pling protocols. This is our goal in the present work: 
here we study bang-bang decoupling in the case where 
the quantum environment can become critical. 

Many-body environments displaying critical behavior 
have been recently investigated in great detail, in order 
to study the sensitivity of decoherence to environmental 
dynamics (see, e.g., [l3|, El])- Close to a quantum crit- 
ical point the environment becomes increasingly slower 
(a phenomenon known as critical slowing down). We 
analyze the decoherence process of a two level system 
(qubit) coupled to an environment modeled as a one- 
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dimensional lattice of spins interacting through an Ising- 
like coupling. We focus on the suppression of qubit de- 
coherence through a bang-bang control procedure, and 
study how the occurrence of a quantum phase transition 
(QPT) in the bath modifies the effectiveness of the con- 
trol procedure. Our analysis is focused on the behavior 
of the Loschmidt Echo (LE) [28] , whose study has given 
new insights into the decoherence process of quantum 
spin chains. We discuss the application of a pulse train 
to the qubit and show its effectiveness in quenching qubit 
dephasing, especially at the critical point. 

The paper is organized as follows. In Sec. [IT] we intro- 
duce the model and pertinent notation. The control pro- 
cedure, based on a sequence of pulses that repeatedly flip 
the state of the system, is described in Sec. lIII| where we 
also derive an explicit expression for the LE in the pres- 
ence of such control. We then provide a detailed analysis 
of its effects in the limiting cases of a single qubit-bath 
link fSubsec. llII A[) and a spin-star model fSubsec. HIIB[) . 
Finally, in Sec.[lV]we discuss our results. In the Appen- 
dices we provide an analytical formula for evaluating the 
LE in the presence of control (App. |A"|) . we perform a 
perturbative analysis in the pulse frequency of the LE 
f App. [B]), and discuss in detail a closed- form formula for 
the LE in the spin-star model (App. IC|) . 

II. MODEL AND NOTATION 

We consider a two level quantum system S (qubit) cou- 
pled to an interacting spin bath E (environment), com- 
prising a linear chain of N spin-1/2 particles, modeled by 
a transverse field Ising model. The Hamiltonian reads 

Ho = Hs + He + Hint j (1) 
where Hs and He are the free Hamiltonians of S and E: 

Hs = -y(l-T*) = -u,o||>UI , (2) 

N 

n * = - J E( ff i CT i+i+ A(7 l); ( 3 ) 

here af and T a (with a = x,y,z) indicate, respectively, 
the Pauli matrices of the ith spin of the chain E and of 
the qubit S, whose ground and excited states are denoted 
by |t) and |j). In this work we will use periodic bound- 
ary conditions, therefore we assume <J% + i = o~f. The 
constants J and A are the interaction strength between 
neighboring spins of the bath and an external transverse 
magnetic field, respectively (in the following, the energy 
and the time scale are taken in units of J, therefore, when 
not specified, we will implicitly assume J = 1). We sup- 
pose that the system is coupled to a given number of bath 
spins [la ]: 

jm 

H int = -e\l){l\®J2 a h ( 4 ) 

3=h 



where e is the coupling constant and m the number of en- 
vironmental spins to which S is coupled. The LE can be 
calculated for a generic sequence {ji, . .. j m } of system- 
bath links. In the following, however, we consider the 
cases m = 1 and m = N. We expect that the generic case 
will be a quantitative interpolation between these two ex- 
tremes but no new qualitative features should emerge. 

With the above choice of Hs and Hint , the populations 
of the ground and excited state of the qubit do not evolve, 
since [t z ,Hq] — 0, and we can study a model of pure 
dephasing. 

As usual, we assume that the initial global state of the 
system is factorized: 

l*(0)> = (o T |T>+c i U»®|G> l (5) 

so that the qubit 5 is in a generic superposition of 
the ground and excited state, while the bath E is in 
its ground state (i.e., \G) is the ground state of the 
Hamiltonian He)- The evolution of such a state under 
the Hamiltonian ([1]) is dictated by the unitary operator 
Uq = e~ lUat and yields, at time t, the state 

|*(t))=c T |T)|^o(t))+c i e-° t ||)|^(t)) , (6) 

where |<p (*)) = er iU ^ \Q) and \<pi(t)) = e~ m ^ \G) are 
the environment states evolved under an "unperturbed" 
and a "perturbed" Hamiltonian given, respectively, by 

Ht=H E , Hi = He + {l\ntot\l) ■ (7) 

The density matrix of the qubit is p — TVe |\&) 
Its diagonal elements are constant, while off-diagonal el- 
ements decay in time as 

p i < [ (t)=p ir (0)e iu ° t D(t), (8) 

with 

D(t) = (<po(jt)\Mt)) = (G\ e'^te-m.i | G) . (9) 

The decoherence of the qubit is then fully characterized 
by the so called Loschmidt echo Co(t) € [0, 1] of the en- 
vironment: 

£<,(*) = |L>(t)| 2 = | (G\ e -i(-H E +{l\n iat \i))t | G ) p _ (1Q) 

The decay of the LE in the model ©-CP with m = N 
(spin-star model) was first studied in detail in Ref. [I3 |: 
an extension to the more general case m ^ TV, and for 
other spin baths - including the XY and Heisenberg 
models - can be found in Ref. [Hj]. It was pointed out 
that the echo decay is enhanced at criticality, due to the 
hypersensitivity to perturbations of the (time-evolved) 
unperturbed ground state [(fait)). Indeed, at criticality 
the perturbation H- m t is very effective at making the un- 
perturbed state l^o(i)) orthogonal to \<pi(t)), thus lead- 
ing to a strong decay of the echo. Away from critical- 
ity, the perturbation is not so effective at orthogonalizing 
|Vo(*)) and \(fi(t)), whence the echo decays more slowly. 
In the following we investigate these effects when a con- 
trol is also present. Details on how to evaluate the LE 
in both the absence and presence of such a control are 
given in Appendix [XJ 
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III. CONTROLLED DYNAMICS 

Quantum dynamical decoupling procedures aimed at 
actively fighting decoherence hinge either on the action of 
frequent interruptions of the evolution or on the effect of 
a strong continuous coupling to an external field. These 
procedures are known to be physically and, to a large 
extent, mathematically equivalent [HI]. Here we focus on 
one possible procedure, based on multipulse control [3l| . 
Let us formally introduce the control scheme as 



H{u ,t) = Ho + H P {uj ,t), 



(11) 



where Tip is an additional time-dependent Hamiltonian 
that causes spin flips of the qubit at regular time intervals 
through a monochromatic alternating magnetic field at 
resonance: 



H P {u ,t) = J2 v{n) (t)\cos(uJo{t~nAt))T x 



+ sin (u)o(t — nAt)) 



(12) 



Here V^ n \t) is constant and equal to V for the en- 
tire duration rp of the nth pulse (i.e., for nAt < t < 
nAt + Tp), At being the time interval between two con- 
secutive pulses. In this work we only deal with n pulses, 
satisfying the condition 2Vrp = ±7T, and suppose that V 
is large enough to yield almost instantaneous spin flips, 
i.e., we take rp <C At. Therefore, in the ideal limit of 
instantaneous kicks of infinite strength (rp — > 0, V — > oo 
such that Vrp — ±7r/2), the effect of each pulse on the 
qubit is simply a flip, that is described by the operator 



Up = ±i t x 



(13) 



The evolution of the initial state ||SJ) under the Hamil- 
tonian (|11|) in one spin-flip cycle [i.e., two flips, from time 
t = to time t\ = 2 (At + rp) ~ 2 At] is dictated by the 
unitary operator 



U c = e~ 2mAt 



and it is such that 



Up Uo(At) Up Uo(At) (14) 



|*(2At)) = -c T e- MoAt |T)e- m i At e-^ At |G) 



(15) 



This is again a pure dephasing phenomenon, so that all 
relevant information is contained in the off-diagonal el- 
ement (jHJ of the system reduced density matrix. The 
behavior of decoherence is then fully captured by the 
LE: 



C(2At) = (G| ( e m i At e m ^ At ) ■ ^-^Atp-iH.A^ | G) 

(16) 

In general, at a certain time t = 2M At + 1, the evolution 
operator of the global system is given by: 



U 



(17) 



Uo(t)[U c ] M if t<At 

^ U (t- At) Up U {At) [U C ] M if t> At 
where M — [^r], [•] denotes the integer part and i = 
t — 2MAt is the residual time after M cycles. It is now 
easy to write down the LE at a generic time t: 



C(t) = 



(G\ ( e 'WiAigi?< T At 



M 



e %H\t . e -iH r t ^ e ~i"Hi At e -iH r At 



\G) 



,iUlAt e m^At 



-i~H.it p-iHf At ( p -iH j. At p -iUf At 



) M \G) 



if t < At 
if t > At 



(18) 



An explicit formula for evaluating the LE, also in the 
presence of pulses, is given in Appendix [Al 

In the limit of short pulse intervals, and when t is an 
integer multiple of the duration of a single spin-flip cycle, 
t = 2MAt, one can show (see Appendix IB) that Eq. (fT5|) 
can be rewritten as 



C(t = 2MAt) 
where 



{G\e ltH ^\G} +MO(At 2 ), (19) 



H ef f = i^[H l ,Hi}=i^[U\Hi n t\l),H E ] (20) 
is an effective Hamiltonian. By noting that (M Hint \i) — 



a ji we h ave : f° r arbitrary A 



N 



3=31 3 = 1 



where 



e e ff 



At 



(21) 



(22) 



This is the renormalized system-bath coupling constant 
in the presence of multipulse control. We notice that H e s 
does not depend on A (which would appear at 0(At 3 ) 
through the double commutator [[H±, Wf], Wf]). There- 
fore, in the small At limit, the criticality of the model 
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can manifest itself only through \G) in the LE expres- 
sion mp . 

In the next two subsections we turn to a numerical 
study of the LE for the cases of a qubit coupled to one 
spin of the chain [m — 1 in Eq. Q], and the spin-star 
model [m = N in Eq. g}]. 



A. Qubit coupled to a single bath spin 

When m = 1, the system-bath Hamiltonian of Eqs. (JTJ> — 
dU) can be rewritten as: 



Ho 



IOUIU- 



ecr 



JV 

J E(« + i+K) ( 23 ) 



and corresponds to a situation in which the qubit is di- 
rectly coupled to only one spin of an Ising chain with pe- 
riodic boundary conditions (the coupled bath-spin qubit 
is assumed for simplicity and with no loss of generality 
to be the first one) . In Fig. [1] we show the behavior of 
the LE in Eq. (JTSJ) as a function of time, for different 
values of the pulse frequency At. The three panels refer 
to different values of the transverse magnetic field A; the 
thick dashed lines represent the case Co(t) with no ex- 
ternal control [At — > oo in Eq. (fl8|) . or simply Eq. IjlOp], 
Here the environment consists of N = 100 Ising spins, 
and the system-bath coupling has been set at e = 0.25. 

We notice a very different behavior as A is varied. 
Away from criticality (i.e., for A = 0.5 Fig. QJa), and 
A = 1.5 Fig. [He)) the LE in absence of control quickly 
reaches its asymptotic (saturation) value Coo, as indi- 
cated by the dashed black lines. Very fast control pulses 
do improve the situation, but only in the sense that this 
asymptotic value becomes slightly closer to unity. In con- 
trast, slow pulses make the situation much worse: when 
JAt is larger than a certain value, the pulses act as an 
additional source of noise and, as a consequence, the co- 
herence decays (exponentially). On the other hand, when 
the chain is critical (A = 1 Fig. [Hb)) and there is no 
control, the LE decays (albeit only logarithmically [HI]), 
as can be seen from the dashed curve. In this case the 
pulses can be very effective, as a control procedure: when 
JAt < 0.375 decay is suppressed. Again, when At ex- 
ceeds this threshold, decay is enhanced. This situation 
is reminiscent of the transition between a quantum Zeno 
and an inverse Zeno effect [40| . 

In Fig.[5]we show the values of the LE at a fixed time t* 
(we performed an average of C(t) for Jt £ [Jt* — 5, Jt*+5] 
in order to eliminate fast oscillations), as a function of 
At. The different curves are obtained for different values 
of the transverse field A. We set Jt* — 25 so that: i) in 
the absence of pulse control and for noncritical A, Co has 
already reached its saturation value Coo ; n) at criticality, 
the minimum of Co(t) is found exactly at Jt* ~ A/4 (in 
this case N — 100) (la ]. In the panel (b), bars denote the 
corresponding value of Co(t*) without external control. 
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FIG. 1: (Color online) Loschmidt Echo as a function of time 
for a qubit coupled to a JV = 100 spin Ising chain, with e = 
0.25. Panels stand for different values of the transverse field: 
(a) A = 0.5, (b) A = 1, (c) A = 1.5; the various curves in 
each panel are for decreasing pulse intervals At, from top to 
bottom. Insets: magnification at small times t (axes units 
are the same as in main panels); notice that, when A = 1, 
frequent pulses suppress decay for JAt < 0.375 (here and in 
the following figures At values are expressed in units of J). 



The behavior at large pulse intervals At is non-trivial 
and rather interesting: we note that the echo has a min- 
imum and has an almost complete recovery, and that 
the LE for A = 0.5 rises higher than for A = 0.9, 1, 1.1. 
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At [J] 



FIG. 2: (Color online) (a) LE as a function of the pulse fre- 
quency At at a given time t*, for different A. (b) Magnifica- 
tion of panel (a) in the highlighted zone; the bars denote the 
corresponding values of Co(t*) without pulsing. Here we set 
Jt* = 25, N = 100, e = 0.25. 




FIG. 3: (Color online) Panels (a)-(b): LE at a fixed time t* 
as a function of the transverse field, for different values of At. 
Panels (c)-(d): rescaled LE, C(t*)/Co(t*). Notice the widely 
different scales in the y axes of (a)-(c) panels (small At), with 
respect to (b)-(d) panels (large At). Here we set Jt* = 25, 
N = 100, e = 0.25. 



The large At regime is non-perturbative (in the sense of 
the perturbation theory of Section [TTT1 and Appendix |B|). 
Nevertheless, the rise of the LE for large At can be under- 
stood as being due to the fact that, after a short transient 
time t, the LE without control saturates around a con- 
stant value (see the black dashed curves in the insets of 
Fig. [TJ or Ref. 15]). Therefore, if the pulse frequency is 
such that At > i, the effect of the bang-bang control pro- 
cedure will be progressively reduced as At grows, until, 
in the limit At —> +oo, it will completely disappear. In 
other words, the detrimental effect of the control for large 
A is offset by the gradual diminishing of its effect as At 
grows, which allows the LE to recover to its saturation 
value. Moreover, as the insets of Fig. [1] show, for A = 1.5 
the saturation is truly at a constant value; for A = 0.5 
the saturation is an oscillation around a constant value; 
at criticality (A = 1) there is a logarithmic decay of the 
LE, but for a finite system size this decay will eventually 
stop and revivals of quantum coherence will appear. The 
oscillation at A = 0.5 explains why this curve rises higher 
than the other curves in Fig.^a); at a time t* w 1.5, the 
uncontrolled LE in Fig. Q] at A = 0.5 is larger than for 
other values of A. 



The panels (a)-(b) of Fig. [3] display £(t*) as a func- 
tion of A, for different values of At. In panels (c)-(d) we 
plot the rescaled quantities, obtained by dividing £(t*) 
by the corresponding value in absence of pulse control, 
£o(t*). The LE has a maximum not at A = 1 but at 
A > 1, while at A = 1 there is an inflexion point. At 
criticality, the rescaled LE displays a cusp. The cusp 
disappears at JAt > 0.375, in agreement with Fig. [IJb), 
where we observed, at the same value of At, an increase of 
the LE when the control is present. A qualitative expla- 
nation of this phenomenon is straightforward: for short 
time pulses, the renormalizcd coupling constant e c g- in 
Eq. ([22]), and therefore the LE, are only weakly depen- 
dent on A at leading order in the perturbative expansion. 
In contrast, the free echo £q has a downward cusp [l5| 
(present also in Fig.^a) for the spin-star case). The ratio 
must therefore display an upward cusp, as seen in Fig. [3] 
Another way to state this explanation is the following. 
For sufficiently small values of At the bang-bang proto- 
col succeeds at effectively eliminating the environment 
action. The only remnant of criticality is then the weak 
signature of an inflexion point seen in Fig.[3Ja). The echo 
of the uncontrolled system, however, is hypersensitive to 
criticality, as indicated by the cusp. On the other hand, 
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when At is too large (Fig. [3^b)-(d)), the bang-bang pro- 
tocol fails at removing the coupling of the qubit to the 
environment, and the controlled and uncontrolled echos 
behave similarly. 

There are other interesting features in Fig. [3j Panels 
(a)-(b) show that the LE rises for sufficiently large A, 
and (c)-(d) show that the ratio between the decoupled 
and free echoes approaches unity for large A. This can be 
understood as being due to the dominance of the uniform 
magnetic field term A J2f=i a j over the transverse Ising 
term J2f=i a j a j+i m Eq- ([23]) . Indeed, in the limit of 
large A, this means that Hi — He [recall Eq. (J7J)], so 
that [H h Hi] ~ and the LE ~ 1 by Eqs. © and (fT0|), 
Thus, at large A, decoupling is not needed to obtain a 
large LE. 

More interesting is the monotonic rise of the LE visible 
in panel (a) as a function of A for J At — 0.1,0.2, in 
contrast to the maximum around A ~ 1.25 for J At — 
0.3. Indeed, panel (c) shows that decoupling makes the 
situation worse for J At = 0.3 and A > 1.25, and a similar 
trend continues in panels (b)-(d). Thus, in our model 
decoupling is fully effective (i.e., for all values of A) for 
J At < 0.2. 



B. Spin-star model 

The "spin-star" model corresponds to the case when 
the qubit is equally coupled to all the spins of the chain 
[m = TV in Eq. Q]. This situation is opposite to the one 
considered in the previous subsection. Interestingly, in 
this limit the model is almost solvable. The system-bath 
Hamiltonian of Eq. |T]) reads: 



LP*) 



N 

Ho = u> ui - •/£ b^+i + ( A + j 



(24) 



We first notice that Hi(X) = H-[(X) = He{X), where 
A = A + e/J. Therefore, both the perturbed and the 
unperturbed Hamiltonians describe an Ising model with 
a uniform transverse field, and can be diagonalized ana- 
lytically by means of a standard Jordan- Wigner-Fourier 
transformation, followed by a Bogoliubov rotation. De- 
tails on how to evaluate the LE of Eq. (fTTj)) for a spin-star 
model can be found in Appendix [Cl where we show that 



(G\e ltH ^\G) = []cos(8te cff A fc ) 



(25) 



fe>0 



with Afc = sin(27rfc/7V) and e c s defined in Eq. ([22]) . In 
the limit of small e c s , while keeping t finite, we can ap- 
proximate this as 

(G\ e ltn ^ \G) ~ J[ e -hi^ a A k f = e -U^)\ (26) 

fc>0 

where we have defined 




T = 64 ^ A 



(27) 



fe>0 



FIG. 4: (Color online) LE for the spin-star model (the param- 
eters of the simulation are Jt* = 10, e = 0.01 and N = 300). 
(a): Dependence on A of the LE without external control at 
fixed time, (b)-(c): LE in presence of pulsed control with 
frequency At. (d): Renormalized controlled LE. 



We notice that the dependence on A in Eq. (f2"6"|) has en- 
tirely disappeared. This means that, to leading order in 
the pulse interval At, dynamical decoupling is not sensi- 
tive to criticality [Eq. flT^} for the LE and Eq. (J25J) lead 
to C » e -r (* e ° ff ) ; we explicitly checked that, for small 
At and at short times, this formula exactly reproduces 
the data obtained from numerical simulations, which are 
completely insensitive to A in that regime]. This is con- 
sistent with the data shown in Fig. [U In panel (a) we see 
the behavior of the LE in absence of control; we notice a 
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slight dip at A = 1. In panels (b)-(c) we study the LE for 
various At; we observe strong similarities with Fig. [3J in 
particular the weak dependence of the LE on A for very 
small At. In panel (d) the rescaled LE again displays a 
cusp. 

It is remarkable how similar the results are for m = 1 
(qubit coupled to a single spin of the chain) and m = N 
(spin-star model). The consistency of these results and 
the analogies between these two opposite situations lead 
us to conclude that general features of the decoherence of 
the qubit under bang-bang control are largely indepen- 
dent of the number of chain spins coupled to it, at least 
when the chain is close to criticality. 



IV. DISCUSSION AND CONCLUSIONS 



are physically equivalent, there are important practical 
differences among them. Future attention will be di- 
rected towards the exploration of these similarities and 
differences in the context of coupling of a qubit to a crit- 
ical system. 
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We have studied the efficacy of pulsed control of a qubit 
when it is coupled to a spin bath. It is well known that, 
without control pulses, the qubit decoheres particularly 
fast in the vicinity of the critical point. The reason for 
this is that the evolution takes the initial state |^(0)), 
in the form of Eq. ((SJ), into a superposition of the type 
IT) l<A)(£)) + II) IviM) an d the two bath states become 
rapidly orthogonal near the critical point. The appli- 
cation of decoupling pulses to the qubit removes the de- 
pendence of decoherence on the criticality of the environ- 
ment. On the other hand, we also found a regime (larger 
interval At between pulses) such that the control can 
increase the effects of decoherence. Away from critical- 
ity the perturbation is not as effective at orthogonalizing 
\(fo(t)) and \tpi(t)), leading to a slow decay of the echo 
and to relatively less effective control. Therefore, we can 
conclude that in general decoupling is relatively more ef- 
fective near the critical point, since there it results in the 
largest enhancement of coherence. 

From the quantum information processing perspective, 
there is another positive message in these results: sup- 
pose we are trying to preserve the coherence of a qubit in 
the presence of a spin bath. Without decoupling we know 
that the spin decoheres particularly fast in the vicinity 
of the critical point. Therefore not knowing whether we 
are close to criticality when trying to operate a quantum 
computer coupled to a spin bath, is a problem. But in 
light of the results presented here, it follows that appli- 
cation of dynamical decoupling pulses removes this con- 
cern: for sufficiently frequent pulses, decoupling works 
independently of the value of the system-bath coupling 
A, so closeness to criticality does not matter. 

Our analytical and numerical calculations suggest that 
these results seem to be largely independent of the details 
of the model of qubit-environment coupling. Indeed, we 
have considered two extreme situations (qubit coupled to 
a single spin of the chain and qubit coupled to all spins in 
the chain), and obtained the same qualitative behavior. 

Finally, a comparison of different control strategies 
(Zeno effect, decoupling pulses and strong continuous 
coupling) [4l| has shown that, although these procedures 
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We explain here how to evaluate the LE for the Hamil- 
tonian in Eq. (JTJ, and then extend some of these results 
to the case of pulsed control, Eq. Ijlip. This technique 
can be easily generalized to the case of an XY spin bath, 
as has been done in Ref. [l5j . 

By means of the Jordan- Wigner transformation 

o+ = c) exp (in 4 C *) ' °*i = 2c h ~ X ' ( A1 ) 
^ fe=i ' 

we first map the Hamiltonians 7i± and H-\ of the spin 
bath onto a free fermion model that can be expressed in 
the form 

H VI = ^C*, (A2) 

where = (c\ . . .c N Ci . . . cn) (cj being the correspond- 
ing spinless fermion operators) and 

C = a z ® A + ia v <g> B (A3) 

is a tridiagonal block matrix with 

A jtk = -J(5k,j+i + <Sj, fe +i) - 2(A + ej)S jt k (A4) 
Bj.k = —J(8k,j+i — Sj,k+i) (A5) 

such that ej = for 7i-\ , while ej = e Sj j j m for . The 
LE can then be evaluated exactly, by rewriting it in terms 
of the determinant of a 2N x 2N matrix (l5j : 

C (t) = |det (l - r + r e tC ^) | , (A6) 

where r is a matrix whose elements r t j = (^1$ •) are the 
two-point correlation functions of the spin chain, evalu- 
ated in the ground state of the Hamiltonian H^. Eq. (|A6|) 
can be obtained from the following trace formula [421 ]: 

Tr ( e r(A) e r (S) ) = det ^ + e A e B) ; (A7) 
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where T(A) = V, A, , a]n , and at, a, are the creation 
and annihilation operators for a fermion particle state i. 

In the presence of pulsed control, in analogy with the 
free evolution case, Eq. pop , we can rewrite the formula 
for the LE in Eq. (fl8|) in terms of the determinant of 
a 2N x 2N matrix. Indeed the trace formula (|A7|) is 
straightforwardly generalized to products of more than 
two operators [42j by using the following identity: 



-in t t -iU 2 t 



-m n t 



det (1 - r + r e- lC ^e- lC2t . 



e - lC "*) , (A8) 



where we supposed that Hk = J2i,j[^^]ij a l a j an d r o = 
r(A/") with Af occupation number operator [i.e. (ro)^ = 

(-00 1 a\ aj | ) ] ■ 
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of short pulses. To simplify the notations, let us define 
A = iH i, B = iHi, and e = At. We consider Eq. (JT5J) at 
integer multiples of a spin-flip cycle, i.e., t = 2MAt: 



C{t) = Tr \\G)(G\(e £A e £B ) M (e 



-eB\M 



(Bl) 



Now recall the (approximate) Lie sum and product for- 
mulas 



Q(e 2 



(B2) 



e £A e eB e- £A e- £B = e^ A ^+Q(e 3 ). (B3) 



Using this we have 



-eA e ~eB\M 



(e eA e ea ) M {e 
= (e eA e eB ) M - 1 [e £2 [ A ' B ]+0( e 3 )](e- EA e- sB ) M - 1 
= (e eA e £B ) M - 2 [e £2 ^e £A e £B +0(e 2 )](e- £A e- eB ) M ' 



r ,eA„eB\M-lf„-eA e -eB\M-l 



+0(e d ){e £A e e ») M - L {e 



(B4) 



We evaluate here the leading order expansion of the LE 
in Eq. (fT8|) in terms of the pulse interval At, in the limit 



Keeping terms only to leading order 0(e 2 ) we can neglect 
the last line, yielding: 



(e £ V B ) M i 



e -eA e -eB )M 



{e eA e e B) M-2 [e e-[A,B] e eA e eB e -eA e -eB 



Q{e 2 )e 



-sA-eB 



-eA -eB\M-2 



eA e eB\M-2r2e 2 [A,B 
eA eB\M-2r 2e 2 [A,B 



1 + 0(e 3 ) + 0(e 2 )(l - s(A + B))]{e~ eA e- eB ) M - 2 



]+0(e 2 )l](e- eA e- £B ) M - 2 



(B5) 



where in the last line we again neglected 0(e 3 ) terms. 
Continuing in this manner we have 



(e eA e eB ) M (e-' A e^ B ) M = e Me ^ + MO(e 2 )l, (B6) 



which yields Eqs. {I9])-(|20 
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The ground state of the Hamiltonian in Eq. (|C1[) is 



\G(X)) 



fc>0 



(y) I00> fc ,_ fc +isin(|)|ll> fe! _ fc 



(C2) 

where 8k = arctan(Afe/efc(A)), and the kets refer to 
fermion occupation numbers in the two modes k and — k. 
Consider now the space 



Here we derive Eq. (|25|) . We first notice that, in the 
spin-star case, both H-\ (A) = He(X) and Hi(X) = He(X) 
can be written in momentum space, by using the Jordan- 
Wigner transformation (jAlj) followed by a Fourier trans- 
form, in this way: 



HeW = 2JJ2 [ek(\)(clc k + cl k c_k) 



k>0 



iA fc (4cL fc - c_ fc c fe )] 



(CI) 



where £fc(A) = A — cos(27rfc/iV) and A& = sin(27rfc/iV), 
and the sum over k runs from 1 to N/2. 



H fe ® H_ fc = Sp{|00) fej _ fe , |01) fci _ fc , |10> fe> _ fe , |H) fci _J . 

Since the subspaces Sp{|00) fe _ fe , | _ k } and 
Sp{|01) fc _ fe , |10) fe _ fe } are not coupled by He, and 
since \G(X)) lives in the former two-dimensional 
subspace, we can rewrite the Hamiltonian over the 
Sp{|00) fe _ k , |H) fc _k\ subspace, up to a constant, as 

HeW = 2jJ2[^(m Z k + = E Ue ^ > 



k>0 



k>0 



(C3) 

where Ej? and "E k generate an SU(2) algebra and are de- 
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fined as 



Now, 





= °i c -k + c -fcC/t, 


(C4) 




= -i(c' k c'_ k ~ C-feCfe), 


(C5) 




= c [ k c k + J_ k c- k - 1 . 


(C6) 



The problem of evaluating (G| 



,WH e {f 



|G) is now reduced 



to computing the action of the 2x2 matrix [H[ , H-\] over 
the subspace Sp{|00) fc _ k , 1 11) fe _ k }. We can rewrite the 
ground state as 



|G(A)) = 



k>0 

®fe>o|G fe (A)) , 



i sin 



(C7) 



where now |0) fc and |l) fe are the standard ±1 eigenstates 
of T, k . Over this subspace, using the fact that Wj.(A) = 



W T (A) + e£ 2 , with £ 2 = Y* =l <?j, we have that 
At 

[ir,'rt T j =4ze c3 

fe>0 



Sg|G fc (A)) = 
so that 



cos 



|l) fe + ^in(y) |0) fe 



<G fc |££|G fe )=0, 



(C9) 



(CIO) 



and 



(G k \e m ^-*\G k ) = (G k \e m ^ A ^\G k ) 
= (G fc |cos(8te off A fc )l 

-isin(8te cfr A fe )I]^|G fc ) 
= cos(8te eff A fc ) (Cll) 



Therefore 



= 8e cff J2 A kX x k = J2 ?W ■ ( C8 ) which is E( l' <E> 



(G\e un ^\G) = Il k>0 (G k \e m ^^ \G k ) 

= n fc>0 cos(8ie cff A fe ), (C12) 



fe>0 



fc>0 
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